COMPRESSIVE INVERSE SCATTERING I. HIGH FREQUENCY SIMO/MIMO 

MEASUREMENTS 



ALBERT C. FANNJIANG 

Abstract. Inverse scattering from discrete targets with the single-input-multiple-output (SIMO), 
multiple-input-singie-output (MISO) or multiple-input-multiple-output (MIMO) measurements is 
analyzed by compressed sensing theory with and without the Born approximation. 

High frequency analysis of (probabilistic) recoverability by the L^-based minimization/regularization 
principles is presented. In the absence of noise, it is shown that the L^-based solution can recover 
exactly the target of sparsity up to the dimension of the data either with the MIMO measurement 
for the Born scattering or with the SIMO/MISO measurement for the exact scattering. The sta- 
bility with respect to noisy data is proved for weak or widely separated scatterers. Reciprocity 
between the SIMO and MISO measurements is analyzed. Finally a coherence bound (and the 
resulting recoverability) is proved for diffraction tomography with high-frequency, few-view and 
limited-angle SIMO/MISO measurements. 



1. Introduction 

A monochromatic wave u propagating in a heterogeneous medium characterized by a variable 
refractive index n = \/l + v is governed by the Helmholtz equation 

(1) An(r)-Fa;2(l-Fz^(r))n(r) =0 

where € C describes the medium inhomogeneities. For simplicity, the wave velocity is assumed 
to be unity and hence the wavenumber uj equals the frequency. 
Consider the plane wave incidence 

(2) n'(r) = e'^""-^ 

where d E S'^~^ , d = 2,3, is the incident direction. The scattered field = u — then satisfies 

(3) Am' + uPv!" = -u?vu 
which can be written as the Lippmann-Schwinger equation: 

(4) u\t) = lo^ f u{r'){u\r') + u%r'))G{r,r')dr' 

where G is the Green function for the operator —(A + u'^) (see Appendix [a| for the reciprocal 
formulation) . 

The scattered field has the far-field asymptotic [35j 

(5) ^^(^) = ^,:^J^(^(^'^) + ^(^))' " = ^ = 

where the scattering amplitude A is determined by the formula |44j 

(6) ^(f,d) = — / iy{r')u{r')e-''^'''-^dr'. 

In the inverse scattering theory, the scattering amplitude is the observable data and the main 
objective then is to reconstruct z^ from the knowledge of the scattering amplitude. In this paper, 
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scatterers 



Figure 1. Far- field imaging geometry 



we use the L^-minimization principle called the Basis Pursuit to study the inverse scattering from 
point scatterers. Note that since u in ([6]) is part of the unknown, the inverse scattering problem 
is nonlinear. Physically speaking, the nonlinearity is the consequence of multiple scattering among 
the different scatterers. 

The standard theory of inverse scattering asserts the injectivity of the mapping from £ 
with a nonnegative imaginary part to the corresponding scattering amplitude for a fixed frequency 
in three (or higher) dimensions (Theorem 5.5 of [13j. See also \2A\ [35l HOl I^H H3] for similar 
results, [291 EH ESI US] for inverse boundary- value problem and [m |30l [311 [35] for inverse obstacle 
scattering ). In this case, the refractive index can be determined uniquely by the full knowledge of 
A(r, d),Vd, f, for a fixed u. Indeed, as A is analytic in both d and r, it suffices to know A for a 
countably many incident and sampling directions in order to determine uniquely. As far as we 
know, the uniqueness in two dimensions with a fixed frequency is still an open question. What is 
known for two dimensions is that the uniqueness holds if the scattering amplitude is given for an 
interval of frequencies [T3| . 

Of obvious theoretical interest, the uniqueness result by itself is of limited practical interest. All 
existing methods for determining the refractive index without linearizing the problem are based 
on the constrained nonlinear optimization in the L^-norm for which the exact recoverability is 
usually difficult to establish, especially in the case of undersampling [13l[Il!- In this paper we show 
that a target is the unique, global minimizer of an optimization principle based on the L^-norm 
if the target satisfies certain sparsity constraint. Moreover, this L^-minimization problem can be 
effectively solved by linear programming as well as various low-complexity greedy algorithms. 

In this paper we focus on the two dimensional setting (r = (x,z) G M^) for the aforementioned 
reason as well as the notational simplicity. Although the details of the results are dimension- 
dependent, our approach is not limited to higher dimensions. We discuss the three dimensional 



case briefly in Sections 2.4 and 2.5 



Consider the medium with point scatterers located in a square lattice £ = {rj = (xj, Zi) : i = 1, m} 
of spacing £. The total number m of grid points in £ is a perfect square. Without loss of generality, 
assume xj = ji£, Zj = where j = (ji — l)^/rn + j2 and ji, j2 = 1, \/m. Let Ujjj = 1, m 
be the strength of the scatterers. Let S = {rj^. = {xi.,Zi^) : j = 1, ...,s} be the locations of the 
scatterers. Hence Vj = 0,Vrj S. When there is no risk of confusion, we shall write v = {vj) in 
the sequel. 

The scattering amplitude for this medium is a finite sum 
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(7) ^(f-,d) = |-^,.,.^(r,)e--^-. 



Moreover, in analogy to (|4j), the exciting field u{ri ) satisfies the Foldy-Lax equation [STj |36] 



(8) u{rii) = u\Yii) + u?^^G{Yii,ri^)vi^u{ri^), / = l,...,s 

where all the multiple scattering effects are included but the self field is excluded to avoid blow-up. 

2. Methods and results 

2.1. MIMO Born scattering. First consider the Born approximation to ([T]). In the Born approx- 
imation (also known as Rayleigh-Gans scattering in optics), the exciting field u{vi.) is replaced by 
the incident field tt'(rj^. ). This approximation linearizes the relation between the the target strength 
and the scattering amplitude and is valid for sufficiently weak or widely separated scatterers. 
For the Born scattering, we define the target vector X = v ^ C™" and use multiple incident waves 

(9) u\{r) = e^^(^^^^6k+xcose^)^ ^ ^ 

where 0^ is the incident angle of the k-th. probe wave. Throughout the paper we consider the 
single-input-multiple-output (SIMO), multiple-input-single-output (MISO) and the multiple-input- 
multiple-output (MIMO) measurements in which for each incident angle 9^ the resulting scattering 
amplitude is measured at the multiple sampling angles Oi,l = l,...,n. After normalization by 
a;^/(47r), the totality of the collected data forms the measurement vector Y G C^". The corre- 
sponding sensing matrix in the linear relationship Y = has the [n[k — 1) + /, j)-entry 



(10) 



-iuj(zj sin di+Xj cos di) \ 



where Oi is the sampling angle of the l-ih sensor. 

Recent breakthrough in compressed sensing has established the insight that the target can be 
recovered exactly with nearly minimum sensing resources by the L^-minimization principle, called 
basis pursuit (BP) 

(11) min||X||i s.t. *X = y 

if the target is sufficiently sparse and the matrix $ satisfies either the incoherence property or the 



restricted isometry property [3l [5l [6l El [U [TT] . The L^-minimization problem (11) can be solved by 
linear programming [21 [HI [11] or by various greedy algorithms [151 EHl HI] ■ 

In this paper we adopt the incoherence approach to analyzing the SIMO/MISO and MIMO 
inverse scattering problems. Previously we have shown in ^21j that suitably designed SIMO and 
MIMO measurements with planar domains (in three dimensions) produce random partial Fourier 
matrix as the sensing matrix which possesses a nearly optimal restricted isometry constant (RIC) 
with respect to the sparsity of the target. Current compressed sensing theory predicts that ( [11[ ) 
yields a superior performance [11|9]. However, with non-planar domains the SIMO/MIMO measure- 
ments seem to produce a rather poor RIC. Hence we adopt the alternative approach of incoherence 
in this paper. The restricted isometry approach will be taken up in Part II for multi-shot SISO 
(single-input-single-output) measurements with non-planar domains for which the framework of 
random partial Fourier matrix can be restored by special sampling schemes. 

To the best of our knowledge the present paper is the first rigorous study of inverse scattering, 
including multiple scattering, in the framework of compressed sensing. Earlier studies |10| [25l (321 
311 nil [53j of related problems largely take the compressed sensing theory for granted and assume 
(explicitly or implicitly) either the incoherence or the restricted isometry property without proof. 

Let us state the perhaps simplest criterion for exact recoverability of the incoherence approach. 
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Proposition 1. [T71 [23] BP reconstructs perfectly any target X of sparsity 
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(12) 



1 

s < - 
- 2 



m(*; 



+ 1 



where the sparsity s = \\X\\o is the number of nonzero components in X and the coherence parameter 
/i($) is defined as 



max 



Proposition [T] implies that the lower the coherence of the sensing matrix is the more massive 
the exactly recoverable target can be. Under the condition (12), a simple greedy algorithm called 
Orthogonal Matching Pursuit (OMP) can provably find the minimizer of ( 11 ) in at most s iterations 
L48J. 

To construct sensing matrices of low coherence let us define the MIMO-sensor ensemble as follows. 
Let the incident angles dh-,^ = l,---,^, be independently and identically distributed according to 
the probability density function E C'^([— vr, vr]); for every incident angle let the sampling angles 
Qi^l = l,...,n, be independently and identically distributed according to the probability density 
function G C'^([— vr, vr]) where /i > is the degree of smoothness. Define supp(/') = {0 : f\9) 7^ 
0} and supp(/'^) = {9 : ^^{0) ^ 0}. We call 9^ £ [— vr, vr] a Blind Spot if there exists a pair r,r' £ C 
such that 



(13) 



(r — r') • (cos 6** , sin 6'* ) I 



In other words, the set of Blind Spots consists of all the angles between the x— axis and r — r', Vr, r' G 
C. 

In Section [3j we prove that the following coherence bound for the sensing matrix with entries 

(To]). 

Theorem 1. Let the sensing matrix $ be given according to the sensor ensemble. Suppose 



(14) 
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Then the sensing matrix pOjl satisfies the coherence bound 
(15) 



Vp J 



X + ' 



with probability greater than (1 — 5)^ where in general x' (resp. x^) satisfies the bound 

(16) x'<Q(l+a;^)-i/2||/||t,oo, 

(17) resp. <ctil + u;ir^/^f%,oo, 

where \\ ■ \\t,oo is the Holder norm of order t > 1/2 and the constant a depends only on t. If, 
however, supp(/^) (resp. supp(/'^)J does not contains any Blind Spot, then x' (resp. x^) satisfies 
the bound 



(18) 
(19) 



resp. < Ch(l + ujiy^yWh^c 



iiriu,. 



■\<h 



d9^ 



where the constant Ch depends only on h. 



Remark 1. Theorem^ along with Proposition^ then imply that any target of sparsity up to 

-1 



(20) 



1 1 

s < - + - 
- 2 2 



V2K 



n 



can he exactly recovered by BP. 

Ifuji is sufficiently large (which is the high frequency limit referred to in the title), the dominant 
term on the right hand side of (20) is 



(21) 



/up 



in view of (16)-(19). 

Note that the high frequency limit for the Helmholtz equation is different from that for the 
Schrddigner equation. The high frequency quantum scattering is essentially linear (without the 
Born approximation) and can be solved by the Radon transform [35]. 



To improve the sparsity constraint (12), Tropp f50] develops an approach in which the recover- 
abihty is only probabilistic in the following ensemble of targets. Let the target ensemble consist of 
target vectors with at most s non-zero entries whose phases are independently uniformly distributed 
in [0, 27r] and whose support indices are independently and randomly selected from the index set 
{1,2,..., m}. 

The following theorem is a reformulation of results due to Tropp ^50j. We refer the reader to 
|22j for the derivation of Proposition [2j 

Proposition 2. Assume the matrix $ has all unit columns. Let X be drawn from the target 
ensemble. Assume that 



(22) 

and that for q > 1 
(23) 



Ai (*)s < 81n 



m 



-1 



re (0,1) 



gins 
21^ 



1/2 



•S 11 ^ 1,9 
+ - * 2 < 

m 



4ei/4- 



Then X is the unique solution of BP with 
norm of $ . 



•ability 1 — 2r — s Here \\^\\2 denotes the spectral 



Remark 2. When the matrix $ has all unit elements, as in (10), the condition (23) becomes 

1/2 

(24) 3 



qhis 



21n- 



+ — Il*ll2 

mp 



< 



1 



4ei/4 



where p = ffrows in 

Proposition [2] calls for the control of the spectral norm of in addition to /u($), in order to 
relax the sparsity constraint from (12) to (22). 

In Section |4] we prove the following spectral norm bound. 

Theorem 2. Under the assumptions of Theorem^ the matrix $ has full rank and its spectral norm 
satisfies the bound 

(25) 

with probability greater than 



I*ll2 < 2m 



(26) 



1-ci 



np — 1 



m 



n{n—l)p{p—l) 



n,p > 2 



for some constant ci > 0. 

In the SIMO case p = 1 we have 



(27) 

with probability larger than 
(28) 



1-ci 



n — 1 



m 



n{n—l) 



n > 2. 



a lower bound for (26) would be 



Remark 3. The probability bounds (26) and (28) are probably far from being optimal. For np ^ m, 

Inp — 1 



1 — cin{n — l)p{p — I] 



m 



which requires m ^ (np)^ to be close to unity. 

When this is violated, we have to rely on Theorem^and Proposition^ which together guarantees 
recovery with probability greater than (1 — 6)'^ but with a higher sensor-to-target ratio. 

Now we are ready to prove the main result for inverse Born scattering. 

Theorem 3. Let the sensors and the target be drawn randomly from the sensor and target en- 
sembles, respectively, and consider the sensing matrix $ of the entries defined by ^-(10). If (I4) 
holds, then the targets of sparsity up to 



(29) 



s< 8 In 



m 



I ■.V2K 

X + 



can be recovered by BP with probability greater than 



n 



(30) 



1-ci 



np — 1 



n(?i— l)p(p— 1) 



ra 



2b \ (1 - 2t - , n > 2, p>2 



for some constant ci > where, for pn ^ s, q can be chosen as 

In m — In e 



72ei/2 Ins ' 



In the SIMO case the probability bound (30) becomes 

n(n— 1) 

(31) ((l-ciJ— ] -25\(1-2t 



n>2. 



m 



Remark 4. Our results can be extended to the case that the sampling angles are not independent 



of the incident angles by adjusting the probability (26) (and hence (30) and ^3^) in the spectral 
norm bound. An important example is the multistatic data matrix with Oj = —Oj,j = 1, ...,n. Such 
a setting is employed in the well known and widely used MIMO imaging scheme called MUSIC 
(standing for MU Itiple-Slgnal-Classification) with p = n \12\ Il6j (see the Conclusion). 

Previously in [22\ , we have applied the compressed sensing methodology to imaging with the multi- 
static data matrix under the paraxial approximation and the assumption that the point scatterers 
lie on a transverse plane. 
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Proof. First, the coherence estimate (15) and the sparsity constraint (29) imphes (22). 



Now the norm bound ( 25 ) imphes ( 23 ) if 



, q\ns\^^'^ 2s 1 



Hence for np ^ s we can choose q in (23) to be 

In m — In r 



72^6 Ins ■ 

Since Theorems [l] and [2] hold simultaneously with probability greater than 

-26, p,n>2 



and since the target ensemble is independent of the sensor ensemble we have the bound (30) for 
the probability of exact recovery. The proof for the SIMO case is the same. 
This completes the proof of Theorem [3| 

□ 

2.2. SIMO/MISO exact inverse scattering. Next we turn to the exact inverse scattering ([T])- 
^ which takes into account all the multiple scattering effects. As the previous simulation shows 
[22j , multiple scattering can severely degrade the performance of the imaging method based on the 
Born approximation. 

Consider both the SIMO measurement with p = 1 and the MISO measurement with n = 1. Note 
that in the MISO measurement, various plane waves are incident upon the scatterers one at a time 
and the corresponding scattering amplitudes are sampled at a fixed direction. By the reciprocity 
of wave propagation in a time-invariant medium, reversing the incident and scattered waves and 
interchanging their roles leave the scattering amplitude unchanged, 

^(f-,d) = A(-d,-f-) 

(see Appendix [A| for a proof). Therefore, the SIMO and MISO cases are equivalent to each other. 

Hereafter we will restrict our attention to the SIMO case. To this end, we will work with the 
alternative definition of the target vector X = {vju{rj)) G C™". The reason for this is that the 
problem then has the appearance of linear system 

(33) Y = ^X 

where the sensing matrix $ has the entries 

(^34^ = ^-iui{zjS\nei+XjCosei) 

and is independent of the incident field. We then apply Theorem |3] with p = 1 and sufficiently 
large n to recover X with high probability. To recover v from X, we observe that as long as 
'u(rj^) 7^ 0, Vj, the support of v is the same as that of X. Indeed, the target strength v = {vj) G C" 
can be recovered exactly by solving the system of nonlinear equations on the target support as 
follows. 

Define the illumination and full field vectors at the locations of the scatterers: 

[/' = (n'(r,J,...,txXr,J)^Ge 
U = (n(r,J,...,tx(r,J)^GC^ 

Let G be the s x s matrix 

G = [{l-5ji)G{vi^,Ti,)] 
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and V the diagonal matrix 

V = diag(z^ii, 
The Foldy-Lax equation ^ can be written as 

(35) U = U' + io^GVU 
from which we obtain 

(36) U = (l-a;2GV)"V' 
and 

(37) X = Vt/ = V(l-w^GV)"V' 

provided that is not an eigenvalue of GV. The exciting field then determines the scattering 
amplitude by ([T]) which yields the (nonlinear) system 

(38) y = *X = *V (I - u^GVy^ U' 
where the sensing matrix entries are given by ( |34[ ). 
Proposition 3. Suppose 

(39) w"^ is not an eigenvalue of the matrix GV 
and 

(40) is not orthogonal to any row vector of (I - uj^GV) . 



X 



Then the solution V of (31) is given by 

where the division is in the entry-wise sense (Hadamard product). In this case, supp(V) = supp(X). 
Proof. Note that 



V (I - co^GV) ^ = (I - oj^VG) ^ V. 



Hence eq. (37) can be written as 

or equivalently 
(42) 



X = (I - cj^vG) ^ VU' 



(I - uj'^VG) X = VU\ 



Solving (42) for the diagonal matrix V entry-by-entry, we obtain (41) which is well-defined if 
(43) cj^GX -|- contains no zero component. 

Since 

u^GX + U' = {I- uj^GVy^ U' = U 

(|43l) follows from diol). 



□ 



Corollary 1. Condition (43) holds and hence (41) is well-defined if 
(44) w^llGVll < 1/2 

where \\ ■ \\ equals the maximum of the absolute row sums of the matrix corresponding to the operator 
norm on L°° . 



Proof. Clearly [/' + u'^GX contains no zero entry if 



(45) w^llG^II < 1 

since every component of [/' has modulus one. As 



|GX|| = ||GV(l-a;2GV) ^ 



( 45 ) follows from ( 44 ) 



□ 



Theorem 4. Suppose p = 1, (14): (39) and hold. Let X he a BP solution for the system (33) 
with the matrix entries (34) according to Theorem^ The formula (41) recovers exactly the target 
of sparsity 



(46) 



s < 81n 



m 



-1 



X 



with probability at least as in (31) and satisfies the bound ply or (19) depending on whether 
supp(/'^) contains a Blind Spot or not. 



Remark 5. The resonance frequency violating (39) is related to the transmission eigenvalue for 
continuous media where an analogous non-resonance condition is also needed to ensure the existence 
and uniqueness of the solution to the inverse scattering proble m [1131 [37] . 

Ifl is an eigenvalue of uj'^GV , the existence of solution for (35) requires that {/' be orthogonal to 
the eigenspace o/w^GV corresponding to 1. Then other physical constraints (such as the minimum 
energy solution) need to be taken into account in order to obtain a unique solution. 

The simplest example for resonance is this: Two point scatterers have the strengths , i^i^ such 
that sign(z/jj) = sign(i/i2) = sign(G*(ri^, rjj). Then 



GV 







12 J 







is real and symmetric and has the positive eigenvalue \J\vi^Vi^G(i:i^^Vi^\. The resonance frequency 
is \vi^Vi.^\~^l^\G(vi^^Xi^\~^l'^ in this case. 



Remark 6. In view of (36), means that U has no zero component. In other words, the target 



is not shadowed by itself in any way. 



ensemble under (39). 



Since the negation of (40) is an algebraic constraint, (40) are satisfied almost surely in the target 



2.3. Stability w.r.t. errors. Here we consider the situation where the measurement or model 
errors are present. In the former case, the data may be contaminated by noise. In the latter case, 
the errors may be due to, for instance, the fact that the targets are slightly off the grid. In this 
case the target vector X does not represent the true targets exactly due to model mismatch and is 
only the best approximation given the model. In either case, the data vector can be written as 



(47) 



Y = ^X + E 



where E represents the errors. In the case of measurement noise, E is independent of the targets 
while in the case of model mismatch, E depends explicitly on the targets. 

Since A(r, d) is an analytic function of both r and d and hence for the given noisy data, in 
general no solution exists to the inverse scattering problem. Even if a solution does exist, it does 
not depend continuously on the measured data in any reasonable norm. 
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To deal with the problem of ill-posedness we consider, instead of the Tikhonov regularization, 
the L^-regularization 

. 1 



(48) 



mm - \\Y 
z 2" 



i + X Z 



where A is the Lagrangian multiplier to be given below. This is the Lagrangian form of the basis 
pursuit denoising (BPDN) [3l [11] and the lasso in the statistics literature [13 [19]. Let X be the 
minimizer of (48). 



The starting point of our analysis is the following result, due to Tropp |49], concerning the error 
bound and the recoverability of the target support in the presence of noise. 

Proposition 4. jJS] Assume that $ has all unit columns and that \\E\\2 < e. 

Suppose /i(*)s < 1/3. Then the minimizer X of (48) with X = 2e is unique and its support is 
contained in supp(X). Moreover, 



(49) 



l^-^lloo < 3+ v/372 e. 



Remark 7. When the matrix $ has all unit entries, as in (34), and the error term satisfies 
\\E\\2 < n^l'^E, then (49) holds with X = 2n£. 

Define the reconstruction of V to be 

X 



(50) 



V = diag 



U' + u}^GX 



in analogy to (41 ). 

Using Proposition |4] we derive an error bound and a sufficient condition under which the support 
of the reconstruction is exactly the same as the original. 



Theorem 5. Suppose \\E\\2 < en^/^ and let X be the solution to (48) with X = 2en. Assume 

(51) ni^)s < 1/3 

and 



(52) 



2 - (3+ y372)e||G| 
Then (50) is well-defined and satisfies the error hound: 

. ^ n -n 2 ( 1 + CJ^ 1 1 Q 1 1 1 1 V 1 1 ) (3 + ^/S/^)^ 

(53) l|V- V|| < ^ ^ ' ' 



< 



1 



a;2(3 + y372)e||G||) 1 
Moreover, supp(V) = supp(X) C supp(X). On the other hand, if 

w^llGVII < l-(!' + v/572).||V-.|| 



l-2co^\\GV\\ 



IGVI 



(54) 



2 - (3+ Y^)e||V-i 



then supp(X) = supp(X). Therefore under (52) and (54), supp(V) = supp(V), i.e. the support of 
the target is perfectly recovered. 

The proof of Theorem [5] is given in Section [5} 



Remark 8. Condition ( [5i[ ) is slightly stronger than (12) and hence the OMP algorithm can he 
used to solve 09] . 



Condition [SW), (54) and (44) clU say in various ways that the scatterers are either weak or far 
apart. 
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Figure 2. SIMO (left) and MISO (right) measurements with single point source 
(left) and single point sensor (right), respectively. 

2.4. Three dimensions. Here we consider the extension of the coherence bound, Theorem [l| to 
the three dimensional setting. The main point here is to demonstrate the decoherence effect due 
to the extra dimension. For simplicity, we will not consider the improved performance as a result 
of avoiding Blind Spots. 

Instead of a square lattice, the computational domain is a cubic lattice of spacing i. Each side 
of the cubic lattice has m^^^ grid points. In three dimensions, the scattering amplitude has the 
same expression except that the sampling direction f = (a, /3, 7) are parametrized by two polar 
angles 6,(j) as 

(55) (5 = cos cos 0, (3 = cos 6 sin (j), 7 = sin0 



Theorem 6. Let the sensing matrix $ be given according to the sensor ensemble. Supp ose 




holds for some constants 6 and K and suppose f\f^ G C^. Then the sensing matrix (10) satisfies 
the coherence bound 

(56) ;,(*)< 

with probability greater than (1 — 5)^ where in general (resp. x") satisfies the bound 

(57) x'<c(l + ^^)-i||/||i,oo [resp. < c(l + ^^)-i ||f ||i,oo) • 

Consequently, the asymptotic behavior ( |21[ ) sets in faster in three dimensions than in two di- 
mension in general. 

2.5. Diffraction tomography: point sensors or sources. Instead of measuring the scattering 
amplitudes, one could measure the scattered field at a set of sampling points and reconstruct the 
targets from the measurement data. Likewise the incident wave may be a spherical wave emitted 
from a point source instead of a plane wave from far field. This is the measurement with point 
sources or sensors. We will focus on the SIMO/MISO settings and discuss both the two and three 
dimensional cases. 

First consider the simple setting with one point sensor at a fixed location measuring the scattered 
fields due to multiple incident plane waves (Figure [2| right) emitted one at a time. This is a MISO 
measurement with one fixed point sensor. By the input-output reciprocity (Appendix [A|) , this is 
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Figure 3. SIMO (left) and MISO (right) measurements with single incident plane 
wave (left) and single sampling direction (right), respectively. 

equivalent to the SIMO measurement of measuring multiple scattering amplitudes due to one point 
source in near field (Figure [2| left). Hence it suffices to analyze the SIMO case. 
Analogous to Q, the scattering amplitude in the direction f is given by 

2 ™ 

(58) A{i;u') = ^^g(r„ro)l^,e-*-^^■^ ^.^r) = G(r, ro) 

where Q is the Green function including the multiple scattering effects, i.e. 



(59) ^(r,ro) = G(r, fq) + ^ z^jt/(rj, ro)G(r, r^), r / r^, = 0, s 

i=i 

(60) ^(rfc,ro) = G(rfc, fq) + ^ ^'j^(rj, ro)G(rfc, r^), = 1, s 

analogous to the Foldy-Lax equation ([s]). Define the target vector X = (Xj) G C"* 

Xj = g{rj,ro)uj, 

the data vector Y = (Y^,) e 

Air 

Yk = -^A{h) 

where = {cos9k,sm9k),k = l,...,n are the incident directions. Then we have Y = ^X where 



the sensing matrix $ = [^kj] is exactly as in (34). 



We can solve this problem exactly as in Section 2.2 by first finding X and then setting 



Xi 



X, 



where we have used (59). In the noisy case (47), we proceed as before and set 



X, 



G(r,-,ro) + c^2 XiG(ro, r,) 

12 



t 



Figure 4. Near- field imaging geometry 



In other words, Theorems |4] and [5] can be immediately generalized to this setting (one point source 
or one point sensor). 

Next let us turn to the more complicated measurement with multiple point sources and one 
sampling direction (MISO, Figure [3| right) or one incident wave and multiple sensors (SIMO, 
Figure [sj left). By the input-output reciprocity (Appendix |A]) , the two settings are equivalent to 
each other. So we focus on the SIMO case below which requires similar but more delicate estimates 
than before. 

We are primarily interested in the imaging set-up called diffraction tomography jl]. First, assume 
the two dimensional setting with the z-axis as the imaging direction (see Figure [s]) . To describe 
the near field measurement we need the Green function for ([T]) 

G{r) = ^Hi,'\u;\r\) 

where Hq^^ is the zeroth order Hankel function of the first kind and admits the Sommerfeld integral 
representation 

(61) H^^\^\r\) = ^ y e^^(l^l^(")+^'^) 

with 

7(a) = 

m 

As shown in Figure [s] the whole space is divided into three domains D U U where the 
infinite slab contains all the scatterers and D~ , are free half spaces. The sensors are placed 
in either D~ (the reflection mode) or (the transmission mode). 

To fixed the idea, suppose the sensors are located in the line segment of length L on {z = 0} C 
symmetrical w.r.t. to the square lattice of the computational domain, see Figure |4j Let Amin and 
Amax) respectively, be the minimum and maximum distances between the line segment and the 
lattice. In two dimensions, 



da 
7(a) 



\/l — a^, \a\ < 1 
iVa'^ — 1, |a| > 1 



(62) A^ax = A/^(i + ^^/^)2 + (A^in + £V^)^ 

13 



Using (61) we can express the scattered field in D U as 

(63) u%r) = f /" — zy(r')n(r')e^'^(^l^-^'l+"(^-^')\ r G U 

Let aj = (^j,0),j = l,...,n where are i.i.d. uniform r.v.s in a square of length L. Let 
Y = (u^{a.j)) G be the data vector. For the SIMO measurement at near field the sensing matrix 
elements, after factoring out w^, are simply 

G(a,-rO = _!|^e*-(7^'+"fe-^')), j = 1, ...,n, I = 1, ...,m. 

For the three dimensional case, we use the Weyl representation formula for the Green function 



where 

7 = i^/a^TW^, + > 1, 

Using (64) we can express the scattered field at z = as 

n^(r) = -—J dr' /^z.(r0^z(rOe-(^^'+-(----')+'^(^-^')), r G Z?" U Z)+. 



Let aj = {ij^rjj^Q), i = l,...,n where {Cj^Vj) i.i.d. uniform r.v.s in a square of length L. 
Assume again that the square aperture is symmetrical w.r.t. the cubic lattice of the computational 
domain as indicated in Figure |4j For this imaging geometry, the maximum distance A^ax between 
the aperture and the lattice is 



(65) A^ax = V ^(^ + ^"^'/')' + + ^^^'^y + (^min + ^^1/3)2. 

Let Y = {u^{aj)) G C"' be the data vector. The SIMO sensing matrix elements , after factoring 
out w^, become 

^ stt^ y 7 

for j = 1, n, / = 1, m. 
Theorem 7. Suppose 

(66) n,<ie^KVrl^ 

where cq depends on the minimum distance Amin between {z = 0} and the lattice (For d = 2, 
ro = 0{- log A^in); for d = 3,ro = 0{A;^J). 
The mutual coherence obeys 

(67) M*) < |G(A^ax)rM^ + 



(68) M*) < |G(A^ax)rM ^ + 4 I ' ^ 



for some constant c (independent of uj > for d = 2 and lo > \ for d = 3), with probability greater 
than (1 — 6)'^ , where A^ax is given by (62) for d = 2 and (65) for d = 3. 
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Remark 9. In view of the presence of the factor |G(Amax)| , we see that Theorem^ is useful 
primarily in the high frequency limit a; ^ 1 with L = 0(1). 

Analogous to Theorems [4] and [5| we have 

Corollary 2. In the absence of noise (e = 0), let X be a BP solution with the sensing matrix 
of diffraction tomography. Under the same assumptions of Theorem^ the formula (41) recovers 
exactly the target of sparsity smaller than 



(69) 
(70) 



1 + |G(A„ 



1 + |G(A„ 



'IjjL 



( V2K ^ c 



d = 2 



d = 3 



for some constant c (independent of uj > for d = 2 and lo > 1 for d = 3) with probability greater 
than (1 - 5f. 

In the presence of noise, Theorem holds. 



3. Proof of Theorem [H Coherence bound 



Proof. Denote i-j = (cos sin 0j), = (cos 6*^, sin 0^)- 
The pairwise coherence has the form 



(71) 



1 

pn 



p n 

jwdfc-(r-r') _ itjfj-(r-r') 



k=l 



where r, r' are two distinct points in the lattice C. Note that the two summations in (71) are of 
the same type. 

Consider the first summation over k = 1, ...,p. Let 



and 



Pk = cos {ujdk • (r - r')), Qk = sin (wd^ • (r - r )) 



Sp = '^Pk, Tp = y^Qfc- 



k=l 



k=l 



Then the summation can be bounded by 

p 



(72) 



k=l 



Juidkir-r') 



< J\Sp - E5p|2 + \Tp - ETpP + JjES-pP + \ETp\^ 



We recall the Hoeffding inequality [26]. 

Proposition 5. Let Pi, Pp be independent random variables. Assume that Pi G [a;, bi],l = 1, ...,p 
almost surely. Then we have 



(73) ^[\Sp- E5p| >pt]<2 exp 

for all positive values oft. 



2+2 



2pH 
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We apply the Hoeffding inequality to both Sp and Tp. To this end, we have bi—ai = 2,\/l = 1, ...,p. 
Set 

t = K/^, K>0. 

Then we obtain 

(74) P [p-^ \Sp - ESp\ > K/^] < 2e-^'/2 

(75) P [p-^ \Tp - ETp\ > K/^] < 2e~^''/^. 

Note that the quantities Sp, Tp depend on r — r' = (xi — xj, Zi — Zj) but they possess the symmetry: 
^^(r — r') = Sp{r' — r),Tp{r — r') = —Tp{r' — r). Furthermore, a moment of reflection reveals that 
thanks to the square symmetry of the lattice there are at most m — 1 different values \Sp\ and \Tp\ 
among the m(m — l)/2 pairs of (r,r'). 

We use (74)-(75) and the union bound to obtain 



maxp ^ \Sp — KSp\ > K/,^ 



maxp-1 \Tp - ETpl > K/y/^ 



< 2(m - 1) • e 

< 2(m - 1) • e-^'/2 

where the factor 4m is due to the structure of square lattice. 
Hence, by fl2^ 

(76) P 



maxp 



.k=l 



giwdfc-(r-r') _ 

k=l 

Similarly we have for the second summation in ( 71 ) 
(77) P 



maxn 



i=i 



< V2K/ 



> (1 - 2(m - l)e^ 



> (1 - 2(m - l)e' 



By (14) the right hand side of (76)-(77) is greater than (1 — Sy 



Consider E 



X]fc=i 6*'^'^''' • The same analysis below applies equally well to E X]?=i 



If P is the uniform distribution over [— vr, vr] or [— 7r/2,7r/2] then 



E 



.fc=i 



,itjdfc-(r-r') 



: pJo(w|r - r I 



where Jq is the zeroth order Bessel function. In general, the exact expression is not available but 
we are concerned only with the asymptotic for (x>|r — r'| S> 1. 
The Bessel function has the asymptotic 



(75 



Xujr) = J {cos iur - -kI A) + Oiiujr) ^)| , > 1. 
V TTwr ^ ^ 



That is, for ^ 1 there exists a constant c > such that 



(79) 

In general, 

(80) 



Jo(u;|r — r'l) < Vr,r'G£, r / r' 



-E 



,fc=i 



,ia;dfc-(r-r') 



27r 



gici;d-(r— r') A 



f{9)d9, d = (cos 0, sine) 
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which is the Herglotz wave function with kernel in two dimensions. By assumption on (80) is 
a finite sum of integrals of the form 

rb 



(81) 



gja;d-(r— r') j^i 



f{9)d9 



with 7^ in (a, b) whose asymptotic for coi ^ 1 can be analyzed by the method of stationary 
phase (Theorem XI. 14 and XI. 15 of [44J). 

Proposition 6. Let gr,r'{0) = d • (r — r')/|r — r'| which is in C°°([— vr, tt]), Vr, r' € C. 
(i) Suppose ^gr,r'{0) ^ 0, V0 G [a, 6], Vr, r' G C. Then for all f G C(}{[a, b]) 



(82) 



Ju\r-r'\g^y{e) A 



no)de 



<CH{l+u\v-v'\)-'^\\r 



h,oo 



for some constant Ch independent of p. Moreover, since {gry : r,r' G C} is a compact subset of 
C^'^^{[a,b]), the constant Ch can be chosen uniformly for all r,r' G C 

(a) Suppose ^5r.r'(^) vanishes at 9^ G {a,b). Since ^5'r,r'(^*) / ^> there exists a constant 
ct,t > 1/2 such that 



(83) 



Mr-r'\gr.r'ie)fi(^0-^d9 



<ct(l + ^|r-r'|)-V2||/i| 



t,oo 



where the constant a is independent of r,r' G C. 
Note that the condition 



d 
d9 



gr,r'{9)^0, 0G[a,6], Vr,r'G£ 



is the same as saying that (a, b) does not contain any Blind Spot. Combining the estimates for 



X 



I 



oja;d-{r-r') 



fi9)d9 



r 



„itjd-(r-r') 



fd9 



using (76)-(77) and the identity 

UV = {U - U){V -V) + tJ{V - ?) + V{U -U) + UV 



we obtain (15) with probability greater than (1 — 6) . 



□ 



4. Proof of Theorem [2j Spectral norm bound 
Proof. For the proof, it suffices to show that the matrix $ satisfies 

(84) ||-***-I„p||2<l 

m 



where I„„ is the np x np identity matrix with the corresponding probability bound (26). By the 



Gershgorin circle theorem, (84) would in turn follow from 

1 



(85) 



i^m< 



np — 1 

since the diagonal elements of $$*/m are unity. 

The pairwise coherence amounts to calculating the expression 



np 
m 



1=1 



1 

m 



1=1 



-iuj{zi sin^+z; cosd) iu){zi sinS'+z; cosd') iuj[zi sinS+z; cosd) —iui(zi sinS'+z; cosd') 
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Summing over ri,l = l,...,m results in finite geometric series in the longitudinal and transverse 
coordinates since they are equally spaced. We obtain 

m 



np 
m 



1=1 



1 

m 



gitj£(cos 9'— cos S+cos 6*— cos 9')y'm/2 ^—iuj{cos8' —cos0+cos0—cos8')y/m,/2 

giuj£{cos0'—cos9+cos9—cos9') 



(86) 



^iu]£{s'm8' —sin9+sin9—sm6)^/m,/2 ^—iuj(s'm0' —s'm9+s'm9—sin9')^/Tn/2 



gi(^^(sin 6'— sin 9+sin 9— sin 9') ^ 

Using the identity |1 — e^'^\ = 2| sin (</>/2)| we then obtain 



(87) 



np 
m 



1=1 



1 


sin 


a;£(cos 9' — cos 9 + cos 9 — cos 9')y/m/2 




m 


sin 


ujI{cos 9' - cos 9 + cos 9 - cos 9') /2 





sin 


ui£{sm9' — sin 9 + sin^ — sm9')y/rn/2 




sin 


uji{sm 9' - sin 6* + sin 9 - sin 9') /2 





Since 9, 9' are independently and identically distributed according to the sine and cosine of these 
variables have the density functions 



9\{t) 



1 



|:/'(arccos t) and 52(*) 



1 



/'(arcsini), 



respectively. Similarly the sine and cosine of 9,9' have the density functions 

g\{t) = ^ f{aiccost) and gl{t) = ^ /"(arcsint), 



respectively. 

Hence the random variables 



Zi = uji{cos 9' -cos 9 + cos 9 - cos 9') /2 G [-2w^, 2ujl] 
Z2 = a;^(sin 9' - sin 6* + sin 9 - sin 9') /2 G [-2ojI, 2ujI\ 
have the density function 

/z. = ^(ff!*5^*5l*<?D(^), ^ = 1,2. 



Since g*g*g*g is bounded in [—4, 4], we have 

Co 



w^>l, 1 = 1,2. 



for some constant cq > 0. 
Define 

(89) 

and note 



C= min min {|Zi — 7rfc|, — 7rA:|} 
9,9' 9,9' fcez 



2C 



sinC>— , Ce(0,7r/2). 
vr 

Hence the probability that {( > b} for small 6 > is larger than 

(1 - ci6)"("-i)p(P-i) 

where the power n(n — l)p{p — 1) accounts for the number of different pairs of random variables 



involved in (89) 
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By the choice 



we deduce that 



np — 1 



m 



1 



1 



with probabihty larger than 



1-ci 



mb'^ np — 1 

n{n—l)p{p—l) 



np — 1 



m 



In the SIMO case p = 1, (89) becomes 
(90) 



^ = minmin {|Zi — vr/cj, \Z2 — irkl} . 



Hence the probabihty that {( > 6} for small 6 > is larger than 

(1 - ci6)"("-i) 

With 



n — 1 



m 



it follows that 



with probability larger than 



t^m< 



l-ci 



1 



1 



np — 1 

-\ n(n-l) 

n — 1 \ 



m 



□ 



5. Proof of Theorem [5j Stability 
Next, we give an estimate for the smallest component of the exciting field vector U = (1 



Proposition 7. // 

(91) 
then 

(92) 



w^llGVll < 1/2 



1 



< 



1 -a;2||GV|| 
1 - 2a;2||GV| 



1 

bo' 



Here for any vector V, V ^ denotes the vector whose entries are the reciprocal of those of V. 
Proof. We write 



with 



(I - w^GV) U' = U'Q{leR) 
R = (C/i)-! (w^GVC/^ + {uj'^GVfU' + ...) 



which converges under (91). Here and © denote the entrywise (Hadamard) product and sum, 
respectively, of two vectors. Hence 



(93) 



1 



[1 - uj'^Gvy^ 



< 



1 



t/i (1 © R) 
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< 



1 



1 - R 



We also have 
(94) 



< (a;2||GV|| +t^4||GV||2 + ...) 



c^^||GV|| 
1 -w2||GVir 



Substituting (94) into (93) we obtain the claimed bound (92). 



□ 



From the Foldy-Lax equation (36) and (92) we have the following lower bound on the exciting 
field vector U 



< l/bo 



(95) 

and hence bo < |u(rj^.)|, Vj. 

Corollary 3. Suppose fJ.{^)s < 1/3 and 

(96) 6o > (3+ \/372) e||V"i. 

Then supp(X) = supp(X). 

Proof. This follows immediately from the fact 



min|X, | = min |fj,u(ri )| > f° > (3 + ^372)6 

j j \\y W 



and Proposition |4j 



□ 



Proposition 8. The vector + oj'^GX contains no zero entry if u}'^\\GX\\ 00 < 1. In particular, 
this is true for the minimizer X of Proposition^^ under the additional assumption 

(97) 60 > ^^(3 + \/372)e||G||. 
In this case, supp(V) = supp(X). 

Proof. The following calculation is straightforward 

(98) \\GX-GX\\oo < ||G||||1 -X||oo < (3+ y372)e||G||. 
Moreover, since 

GX = GV{I-u}^GVy^U' 

we have the estimate 

IIGVII 



(99) 
Hence 



1 -w2||Gvir 



(100) 



w^||GX||oo < a;^||GX||oo + u;2||GX-GX||oo 
= l-bo + u}^i3 + ^/3/2)e\\G\\ <1 



under the additional condition (97). 

The proof of Theorem [5] can now be completed as follows. 
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□ 



Proof. First of all, (52) is equivalent to (97) and by Proposition [s] the formula (50) is well-defined. 
Subtracting (41) from (50) we can estimate as follows: 



IV- VI 



U'QiX-X)+ Lo^iX -X)QGX - oj^X G{X - X) 



< 



{U' + uj^GX){U' + u^GX) 
1 + uj^\\GX\\^) \\X - X||oo +a;2||X||oo||GX - GX\ 



1 -w2||GX||oo 

where we have used the identity 

(101) U = io^GX + U\ 



X IIC/^^I 



By (|9]J) and (|99|) we find that 
(102) 



u^WGXW^ < 1 



And, since X = (i/jU^rj)), 

\\X\\oo < ||V||||C/||oo. 

This, (95), (98) and Proposition |4] lead to the bound 

(2 + a;2||G||||V||||;7||oo)(3 + 7372)6 



(103) 



V- V < 



60(60 -w2(3 + Y^)e||G 



In view of (101) and (99) we have the following bound 
(104) \\U\\oo <l+uj^\\GX\\^ < 



1 -w^llGVl 



< 2. 



The claimed result (53) now follows from ([iM]) and ([l04]). 

Since (54) is equivalent with (96) it follows from Corollaryjsj Propositionsjsjandjsjthat supp(V) 
supp(V). 

6. Proof of Theorem [6] 
Let d = (a, /?, 7) be parameterized by the angles 6, cp as 
(105) a = cos 9 cos (j), f3 = cos 6 sin 0, 7 = sin 6. 

The pairwise coherence has the form 



□ 



(106) 



1 

pn 



k=l 



("fe,/3fe,7fc)-(r-r') ^ gjaj(aj,/3j,7j)-(r-r') 



where {ctk, f3k,^k), k = l,---,p and {aj, f3j,^j),j = l,...,n are independently and identically dis- 
tributed in the unit sphere according to f^{9,(j)) and f^{6,(j)), respectively. 

The main difference between two and three dimensions is in evaluating the expectation of 

ELi e''^^"*'^''''^''^'^'""'"''' and n"^ Yl]=i e''^'''^"^"'''^'^'''''"' which amounts to calculating the in- 
tegrals 

f7r/2 

d9f\{6 + ^0) cos 9 exp [ia;|r — r'| sin 0] 

-7r/2 
I-tt/2 

/ d9fl{9 + ^o) cos 9 exp [w|r - r | sin 9] 
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(107) 
(108) 



for some Oq depending on r — r' where /| and ff are the marginal density functions 

d<Pf{0,4>) 



fm 



dcf>ne,<p) 



If /} = ff = l/vr, the integrals (107) and (108) become 



2 sin (a;|r — r'| 



UJ\r 



0(1), .«»L 



For the general case, integrating by parts with (|107|) and (|108|) produces 
i 



(109) 
(110) 



uj\r — r 

i 



uj\r ■ 



/i(e + 0o)e^'^l'-'''l''''^ 
/f(0 + 0o)e*^''""'"''''"' 



7r/2 


1 g'ta;|r— r'l sin6 


d 




/-7r/2 


de 


7r/2 


i-n/2 

1 gia;|r— r'l sin 




-77/2 


J-it/2 





from which we obtain the bound 

f.7r/2 



(111) 

(112) 



d0/j(^ + ^o)cos^exp [ia;|r-r'|sin^] 

/■TT/Z 

/ (i6'/^(6' + 6*0) cos6'exp [i(u|r - r'l sin^ 

J -it/2 



-■it/2 
it/2 



< 



f[{9 + 9o)de 

fl{e + eo)dd 

ii/iii,oo 



< 



1,00 • 



7. Proof of Theorem [7] 
7.1. Two dimensional case. The pairwise coherence is a ratio of the numerator 



(113) 

and the denominator 
(114) 



j;G*(a,-r)G(a,-r') 

1/2 



1/2 



'm2 



i=i 

with r = (x,z),r' = {x',z') being two distinct elements of {r^,/ = l,...,m}. Note that G*{aj — 
r)G{aj — r') are i.i.d. random variables since aj are. 

To analyze the summation in (113) we follow the argument of the proof of Theorem [Tj Define 
the random variables Pj,Qj,j = 1, n to be the real and imaginary parts, respectively, of G*{aj — 
r)G{aj — r'). We apply the Hoeffding inequality to both Sn = Z^"=i Pj a-^id T„ = Yl]=i Qj- 

Note that the ranges tq = bj — aj of Pj and Qj depends on the minimum distance Amin between 
{z = 0} and the square lattice (Figure|4]). Indeed for small distance Amin <C 1, ro = 0(— log Amin)- 
We have 



(115) 



1 

max — 

i^j n 



>(i - 5y 



^ [G*(a,- - r)G(a,- - r') - E (G*(a, - r)G(a, - r'))] 



if (pi) holds. 
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Using (61) we have the fohowing calculation for the expectation: 



E{G*(a,-r)G(a,-r')} 
167r2 J -f* y I 



(116) 



167r2 



da da' 



M-zY+z'-y')^iiu(ax-a'x') ^ ^in {{a' - a)ujL/2) 

uLia' — a) 



We divide the integration (116) into three parts: = 1,2,3, the integrals over, respectively, 

Ri = {|a| > \/2}U{|a'| > V2}, 
i?2 = R\r\{\a-a'\<H/{ojL)], 
= Rlr\{\a-a\> H/{u:L)] 

where H <^ ujL will be determined later. 
To estimate /i, note that 

2sin ((a' — a)/e) ^, , n . , ^ i- -i • 

lim ; = d(a — a) m the sense oi distribution. 

£^0 a' — a 

Consequently ujL ■ \Ii\ tends to 



1 



167r2 



da 

\a\>V2 ItP 



^-uj\z+z'\\-y\ ^iuja{x-x') 



< 



-2ajAniin 

167r2 



da 
V2W 



ojL — )• 00, 



and hence is bounded uniformly in cjAmin > 0. 
For I2 and I3 we have the rough estimates 



\h\ < 
\h\ < 



1 



167r2 
1 

167r2 



dada' 



O 



R2 m\ 

dada' 1 



H 

(jjL 



Ra 



o 



where the singularities are inverse-square-root and hence integrable. By choosing H 
obtain 



'LoL we 



(117) 



|E{G*(aj -r)G(aj -r')}| < 



ci 



'ujL 



for some constant ci > 0. 



(118) 



A simple lower bound for the expression (114) is given by 

2 



min |G(a — r) 



where a runs through the entire square aperture of length L and r runs through the entire square 
domain of length i^/m. Since |G| is a decreasing function of the absolute value of its argument, ( 118 ) 
is achieved by the pair of points furthest away from each other. In the assumed imaging geometry. 



Figure [4| the largest distance A^ax between any pair of points is given by (62). Combining (115 
(117) and (118) we obtain that (67) as claimed. 
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7.2. Three dimensional case. We now need to estimate (113) and (114) with r = {x,y,z),r' = 
{x' , y' , z') being two distinct elements of {r^, I = 1, ...,m}. We proceed as before. 

Define the random variables Pj,Qj,j = 1, n to be the real and imaginary parts, respectively, 
of G*{a.j — r)G{aLj — r'). In three dimensions, the ranges tq of Pj and Qj are 0{A'^^^). 

Using (64) we have the following calculation: 



^ E{G*(aj -r)G(aj -r')} 



1 f dadf3 da' dp' 



64tt* J 7* y 
1 f dadp da' dp' 



jg|g-ja;(z7* + (Cj-a;)a+(»,,-2;)/3)gia;(2'7'+(?j-x'')Q'+(»?,-y')/3')| 



^iu}{-Z'^*+z''y') ^iu}{ax-a'x') ^iu}{l3y-(5'y') 



647r4 J 7* y 
2 sin {{a' - a)ojL/2) 2 sin {{P' - P)ujL/2) 
ujL{a' - a) ujL{P' - P) ' 

We divide the integration into several parts as follows. Denote 

Ai = {|a| > 1}U{|q'| > 1}, 

A2 = {\a-a'\> H/{loL)}, 

Bi = {|/3| > 1}U{|/3'| > 1}, 

B2 = {\P-P'\>H/iioL)}. 

Let Ii,i = 1,2,3,5,6,7, be the integrals over, respectively, Ri = Ai Ci Bi,R2 = A^ Ci Bi,R3 = 
Ai n Bf, R4 = A'in B1. 

11 can be estimated straightforwardly by 

^ - 647r4w22,2 |^|2 

for a constant c and is bounded uniformly in w > 1. 

12 can be further decomposed into two parts parts /2j,i = 1,2 corresponding, respectively, to 
integration over R21 = n i?i n A2, R22 = A\r\Bir\ A2 and estimated by 

c [ dadpda'dp' (I H\ 1 

h - /2i + /22<^y |^^,| [h + z;z)-^ 

for a constant c. 

/s can be further decomposed into two parts Is?,, i = 1,2, corresponding, respectively, to integra- 
tion over R^i = ^1 n 5f n B^, R32 = AinBfnB2 and estimated by 

c f dadpda'dp' (I H\ 1 

h = /3i + /32<^y |^^,| + 

for a constant c. 

Ii can be decomposed into four parts 1^,1 = 1, 2, 3, 4 corresponding, respectively, to integration 
over 

i?4i = A\r\A2r\Blr\B2, 
i?42 = A\r\A2r\Blr\B^, 
R43 = AlnA^2^Bfr\B2, 
R44 = AlnA^nBlnB^, 
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and estimated by 



< 



647r4 



dadl3da'dl3' /II I H HI H H 



for some constant c. 

Setting H = (cjL)^/^ we obtain 



for some constant C2 > 0. 

Using (118) for the lower bound of (114) we conclude (68) in Theorem [tI 



8. Conclusion 

We have analyzed the SIMO/MISO and MIMO inverse scattering problems by compressed sens- 
ing theory to shed new light on this problem with distinguished history \13\ [BUI [5^ . We have 
obtained several main results: Theorem [s] concerns the recoverability by minimization with the 
MIMO measurement under the Born approximation, Theorem |4] addresses the recoverability by 
minimization with the SIMO/MISO measurement including multiple scattering and Theorem 
[5] asserts stability to measurement or model errors for weak or widely separated scatterers under 
a stronger sparsity constraint. We have also analyzed the diffraction tomography with few views 
and limited angles and proved the coherence bound (Theorem [7|) analogous to Theorem |3j 

A main limitation to our approach is the assumption of finite (albeit high) dimensional targets. 
Also, the reconstruction succeeds only probabilistically. These limitations are intrinsic to the 
current formulation of the compressed sensing theory which is still evolving and in this regard the 
present paper is only a first step in the new direction. 

On the other hand, the compressed sensing approach is constructive and treats the uniqueness 
and the reconstruction in a unified way. Indeed, the main advantage of this approach is an explicit 
and efficient method (i.e. the Basis Pursuit with the SIMO/MIMO or MIMO sensing matrix) for 
reconstructing the scatterers from the scattering amplitude. Moverover, the aperture can be rather 
arbitrary and the dimension of measurement can be as low as comparable to the target sparsity 
(up to a log (m)-factor). 

It may be worthwhile to compare our imaging method with the MUSIC algorithm which employs 
multiple sensors to collect the n x n multistatic response data matrix where n is the number of 
transmitters/receivers [12| 146 j. When the measurement is carried out in the far field, the 
entry of the response matrix is the measured scattering amplitude for the sampling direction / 
and the incident direction j. It is not known if MUSIC can recover the target support exactly for 
nonlinear inverse scattering. Only the case for the Born approximation has been shown capable 
of exact recovery of the target support in the absence of noise [.27] (see the corrected argument 
in Theorem 4.1, |28J). And the estimate for the required dimension of the measurement for the 
exact recovery is hardly optimal. This result should be compared to Theorem [3] with p = n and 



6j = —Oj,j = 1, ...,n, in particular the sparsity constraint (29) for compressed sensing versus the 
necessary condition n > s for MUSIC. This represents a significant reduction in the number of 
sensors when the sparsity of the target vector is large. 

In a separate paper [20] we propose novel multi-shot single-input-single-output (SISO) com- 
pressive imaging methods and demonstrate their superior performances including the capability of 
imaging extended targets. We also present in [20] numerical comparative study of the respective 
performances of the SIMO/MIMO and multi-shot SISO schemes. 
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Appendix A. Input-output reciprocity 
More generally, consider the Helmholtz equation with a source 

Au(r) + u?{l + v{y))u{y) = -/(r) 
which can be solved by using the Green function Q as 

^x(r) = j g{vy)f{v')dv'. 

By slight abuse of notation, we shall write the solution as 

u = gf 

where Q stands also for the corresponding propagator. 
Because the incident wave is governed by 

AuXr)+a;V(r) = -/(r) 

we can write 

u{v) = -g [Au' + (r) 

= g[-{l^ + u?{l + v))u' + u^vu']{v) 
= u\r)+u^g[vu']{r). 
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Hence the scattered wave = u — can be expressed as 
(119) u%r)=uj^ [ g{r,r')u{r')u\r')dr' 



which is the reciprocal representation to Q. 
For point scatterers, (119) becomes 



(120) 



u 



(r) = ^ i^j^(r, rj)M'(rj), r / r^, k = 1, s. 
Substituting the Foldy-Lax equation 

(121) g{r,rj) = G{r,rj) + u)^'^UkG{rj,rk)G{r,rk), r ^ rj, j = I, s 



in (120) we obtain 

m m 

it^(r) = uj^^Vju'{rj)G{r,rj) +uj'^^^UjUkg{rj,rk)u'{rj)G{r,rj] 

j=l j=l kj^j 

which can be rewritten as 



u 



(r) = 5j^k^^'^Uku\Yk)G{r, rj) + X rfc)^'(rfc)G'('^> ^^j. 

j=l k=l j=l kj^j 

m m 

(122) = X^X^ [(5j,fca;^i/fe + (1 - 5j,fc)a;^z/jffc^(rj, rfc)] n'(rfe)G(r, r^) 

j=i k=i 

(see [331 El] for a similar, but slightly erroneous, expression). 



To obtain the alternative expression for the scattering amplitude, let r — )• oo in ( 122 ) and extract 
the plane wave spectrum by (61) for d = 2 or (64) for d = 3. The scattering amplitude in the 
direction r is given by 

^ m m 

(123) Aii-, u') = —^Yl [hkoo^^k + (1 - 5,,k)oo\ukQ{rj,rk)\ u\vk)e-'^'-'^ 



j=i k=i 



where v} is not necessarily a plane wave. 

In the case of a plane wave incidence ^ we observe the symmetry between the incident and 
scattered plane waves in (123). Therefore, reversing and interchanging roles of the incident and 
scattered waves do not affect the scattering amplitude, i.e. ^(f,d) = yl(— d, — r). This is the 



reciprocity referred to in Section 2.2 



In the case of a point sensor located at fq and an incident plane wave, the measurement data is 
given by 

m m 

(124) n^(ro) = [^i.k^^^k + (1 - 5j,k)i^%Vkg{vj,Vk)\ e^"^-'-^G(ro, r,) 

j=i k=i 

whose right hand side can also be interpreted as the scattering amplitude in the direction — d when 
a point source is placed at ro, i.e. A(— d, n') with u^{r) = G(r,ro). This is the reciprocity referred 
to in Section [2?5l 
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